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Abstract 

Changes  in  climatic  parameters  are  important  in  that  they  affect  the 
military’s  ability  to  perform  its  national  defense  mission  and  to  manage 
military  lands.  The  military  must  have  the  ability  to  predict  climatic 
changes  on  specific  installations.  To  support  this  need,  the  research  tested 
and  evaluated  the  application  of  six  multivariate  approach  techniques  to 
predict  climatic  changes  on  a  specific  Army  installation,  Fort  Benning,  GA. 
The  six  approaches  were  tested  for  their  ability  to  identify  where 
anticipated  future  conditions  might  be  found  today  using  a  set  of  19 
bioclimatic  parameters  derived  from  climate  change  data.  The  evaluation 
found  that  the  Primarily  Analogous  Multivariate  approach  developed 
during  this  research  clearly  distinguished  itself  from  the  other  five 
approaches  in  that  it  successfully  determined  future  climatic  factors  at  the 
installation,  as  well  as  at  the  sub-installation  level  using  the  currently 
available  climatic  data. 
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1  Introduction 

1.1  Background 

Climate  model  projections  summarized  in  the  2007  Intergovernmental 
Panel  on  Climate  Change  (IPCC)  report  (IPCC  2007a)  indicate  that  global 
surface  temperature  is  likely  to  rise  between  1.1  and  6.4  °C  during  the  21st 
century.  In  February  2010,  in  response  to  climate  change  forecasts,  the 
President’s  Council  on  Environmental  Quality  (CEQ)  issued  draft  guidance 
to  all  Federal  agencies  concerning  the  manner  in  which  climate  change 
should  be  included  in  the  evaluation  of  environmental  effects  under  the 
National  Environmental  Policy  Act  (Sutley  2010).  Specifically,  the 
guidance  states  that: 

Agencies  should  consider  the  specific  effects  of  the  proposed  action 
(including  the  proposed  action’s  effect  on  the  vulnerability  of  affected 
ecosystems)... 

The  Quadrennial  Defense  Review  (QDR)  was  the  first  high-level 
Department  of  Defense  (DoD)  publication  to  address  the  issue  of  the 
growing  need  to  consider  risks  and  response  strategies  for  climate  change 
(DoD  2010).  In  the  QDR  the  DoD  explicitly  acknowledged  that  climate 
change  will  likely  affect  the  nature  and  scope  of  future  missions,  as  well  as 
training  and  testing  assets  of  military  installations. 

Installation-level  Integrated  Natural  Resources  Management  Plans 
(INRMPs)  provide  an  appropriate  context  for  considering  such  climate 
change  impacts  (Legacy  2009).  Natural  communities  require  decades  to 
mature  and  evolve.  At  the  installation  level,  the  purpose  of  the  INRMP  is 
to  plan,  and  to  coordinate  with  other  stakeholders  in  the  region,  landscape 
or  ecosystem  level  management  over  long  time  periods.  The  goal  of 
INRMP  is  to  ensure  both  good  stewardship  and  no  net  loss  of  the  lands 
that  support  the  military  mission. 

To  ensure  that  the  Army  will  meet  its  training  and  materials  testing 
requirements  in  the  future,  it  must  be  able  to  anticipate  climate-induced 
changes  that  may  affect  the  activities  at  each  of  its  major  installations.  To 
do  this,  the  most  basic  need  is  to  identify  how  the  environmental 
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conditions  that  support  the  local  ecosystems  might  shift,  especially  if  that 
shift  favors  a  different  (i.e.,  changed)  ecosystem.  However,  most  literature 
examining  the  effects  of  climate  change  on  military  installations  has  dealt 
with  the  effect  of  rising  sea  levels  on  coastal  areas,  a  concern  more 
important  to  Navy  and  Marine  interests  than  to  Army  installations. 
Installation  land  managers  and  planners  need  scientifically  reliable 
information  that  addresses  climate  change  in  the  context  of  the 
ecosystems  in  which  they  reside. 

To  meet  this  need,  the  Army  formed  a  climate  change  working  group  to 
undertake  a  high-level  climatic  change  assessment  to  better  target  military 
concerns.  *  The  initiative  is  intended  to  assess  vulnerabilities  across 
locations  and  missions  at  a  regional  scale  and  to  integrate  climate  change 
considerations  into  installation  plans  and  assess  technical  assistance 
needs.  Specifically,  the  effort  will  help  land  managers  maintain  installation 
lands  if  (and  when)  the  ecosystem  changes,  while  still  supporting  the 
installations’  military  mission.  This  study  is  one  in  a  series  (Lozar  2011, 
2012a,  2012b;  Westervelt  2011)  that  focuses  on  those  needs  at  individual 
military  installations.  This  work  analyzes  forecasted  ecosystem  shifts 
around  Fort  Benning,  GA  and  attempts  to  predict  the  geographic  locations 
of  combinations  of  environmental  conditions  and  the  ecosystems 
associated  with  those  conditions. 

1.2  Objective 

The  overall  objective  of  this  project  is  to  explore  approaches  that  can  help 
installation  natural  resource  managers  anticipate  ecosystem  stresses  and 
shifts  that  may  occur  in  response  to  continued  climate  change.  The  specific 
objective  of  this  work  was  to  identify  geographical  areas  that  currently 
experience  environmental  conditions  that  are  the  same  as  climate 
conditions  forecasted  to  occur  at  target  location  in  the  future.  A  secondary, 
related  objective  was  to  determine  whether  currently  available  climatic  data 
were  sufficient  to  successfully  analyze  climate  change  at  the  sub-installation 
level  with  the  currently  available  climatic  data. 


*  http://www.army.mil/article/101844/Climate_Change_Task_Force/ 
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1.1  Approach 

The  objectives  of  this  work  were  met  in  the  following  steps 

1.  Climatic  data  was  obtained  from  various  authoriative  sources.  These 
consisted  of  19  bioclimatic  parameters  that  represented  Global  Climatic 
Models  (GCMs)  and  three  future  scenarios  for  each  of  the  time  periods 
centered  around  the  dates  1990,  2025,  2055,  and  2085.  A  consensus 
database  was  created  from  these  sources.  This  database  formed  a  solid 
foundation  that  made  it  possible  to  reliably  track  the  predicted  changes. 

2.  Four  Multivariate  Analysis  (MVA)  techniques  were  applied  to  climate 
change  data  pertaining  to  the  area  around  Fort  Benning,  GA  to  determine 
the  relative  ease  and  efficacy  of  the  different  techniques. 

3.  The  results  from  the  application  of  the  techniques  were  compared  and 
evaluated. 

4.  Conclusions  were  drawn  and  recommendations  made  for  future  research. 

1.3  Scope 

This  investigation  inspects  the  data  and  analysis  techniques  already 
available.  No  attempt  was  made  to  generate  new  climate  change  preditions 
techniques  or  to  program  new  software.  This  work  used  scientifically 
validated  climate  change  data,  not  weather  or  weather  extremes,  and 
assumed  that  the  military  missions  at  installations  will  remain  the  same  as 
they  are  today. 

1.4  Mode  of  technology  transfer 

It  is  anticipated  that  the  results  of  this  work  will  be  used  to  further 
development  of  techniques  to  predict  future  climatic  factors  on  the 
installation  and  sub-installation  level. 
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2  Description  of  Climate  Change  Data  Used 

2.1  General  concerns  in  the  choice  of  the  dataset  adopted 

Research  scientists  are  particularly  sensitive  to  the  issue  of  choosing  data 
biased  toward  a  particular  predisposition.  Although  military  groups  do  have 
“favored”  datasets  for  climatic  change  research,  this  research  took  great  care 
to  follow  guidelines  in  choosing  data  that  provided  the  greatest  viability  and 
overall  acceptance.  Those  guidelines  include  using  data  that  were: 

•  generated  by  others  (outside  the  local  research  group)  with  specialized 
expertise 

•  generated  in  a  documented,  scientifically  accepted  manner 

•  freely  available  to  ensure  easily  and  accurately  repeatable  research. 

These  guidelines  were  adopted  not  only  to  ensure  that  the  results  of  this 
work  will  be  scientifically  valid,  but  also  to  ensure  that  they  will  be 
compatible  with  any  dataset  eventually  the  military  chooses  to  adopt. 

2.2  The  Intergovernmental  Panel  on  Climate  Change 

The  most  respected  GCMs  generate  predictions  based  on  a  set  of 
conventions  disseminated  through  the  IPCC.  Such  standardization  is  meant 
to  facilitate  comparison  between  models.  The  IPCC  reports  (2007b)  are 
intended  to  reflect  the  scientific  consensus  among  the  experts  in  the  field. 

All  data  used  in  this  research  follow  the  IPCC  standards. 

2.3  The  climate  change  scenarios 

The  IPCC  has  established  a  series  of  standard  future  scenarios  to  assist 
with  coordination  and  comparison  between  modeling  efforts.  This 
international  standard  set  of  scenario  types,  named  the  “SRES”  (after  “The 
Special  Report  on  Emissions  Scenarios”),  was  prepared  by  the  IPCC  for  the 
Third  Assessment  Report  (TAR)  in  2001  on  future  emission  scenarios  to 
be  used  for  driving  GCMs  to  develop  climate  change  scenarios.  The  SRES 
were  also  used  for  the  Fourth  Assessment  Report  (AR4)  in  2007. 
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This  work  used  these  SRES  emission  scenarios: 

•  Ai(B):  globally  homogenous  rapid  economic  growth 

(with  B  variation  =  a  balanced  usage  of  both  fossil  and  non¬ 
fossil  fuel  energy  sources) 

•  A2:  locally  heterogeneous,  regionally  oriented  economic  growth 

•  Bi:  globally  homogenous  sustainable  economic  growth. 

2.4  The  major  climate  models 

About  21  major  GCMs  were  used  in  the  development  of  AR4.  This  research 
chose  to  use  the  models  that  have  had  the  greatest  number  of  validation 
studies  and  the  longest  period  of  development  (1  to  2  decades): 

1.  GFDL  Model  -  NOAA  Princeton  (gfdl_cm2_i) 

2.  NASA  GISS  (giss_model_er) 

3.  United  Kingdom  Hadley  Model  (ukmo_hadcm3) 

4.  Canadian  (CCCma)  Model  (cccma_cgcm3_i__t47) 

5.  NCAR  Boulder  (ncar_ccsm3_o) 

6.  Australian  Model  (csiro  mk3  5). 

Note  that  the  Fifth  Assessment  Report  (the  AR5)  is  in  progress  at  the  time 
of  this  writing.  While  the  data  are  just  becoming  available,  the  types 
available  are  not  nearly  as  mature  as  those  based  on  the  AR4  work  (IPCC 
2013).  For  this  reason,  the  most  recent  model  results  (the  AR4)  were  used  in 
this  study. 

2.5  Characteristics  of  the  data 

The  GCM  data  used  here  have  been  downscaled  or  refined  from  their  initial 
resolution  by  integrating  the  data  with  more  local  concerns  such  as  topog¬ 
raphy,  surface  winds,  evaporation,  and  local  precipitation  (Archive  Collabo¬ 
rators  2013).  Downscaling  using  statistical  approaches  produces  maps  with 
patch  resolution  down  to  30  arc-seconds  (~8oo  meters)  on  an  edge.  Data  at 
this  resolution  were  used  in  the  analyses  described  in  this  report. 

The  data  used  were  averaged  over  a  30-year  period.  For  simplicity,  this  report 
refers  to  the  entire  period  using  the  midpoint  date  of  the  period.  For  example, 
the  period  of  2010  to  2039  will  be  referred  to  as  “2025.” 
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2.6  “Current”  and  “predictive”  data 

The  WorldClim  dataset  (Hijmans  et  al.  2005),  which  represents 
downscaled  data  from  weather  stations  averaged  over  a  period  of  1950- 
2000,  was  used  to  represent  “current”  conditions  (WorldClim  2013).  This 
dataset  contains  global  maps  of  the  minimum,  mean,  and  maximum 
temperatures  and  rainfall  by  month  over  the  period  of  1950-2000.  It  is 
based  on  interpolation  of  weather  station  data  from  across  the  globe  that 
have  been  adjusted  for  the  effects  of  elevation.  This  report  uses  the  term 
“current”  to  refer  to  this  dataset. 

For  the  IPCC  analyses  the  various  GCMs  are  run  using  inputs  defined  by 
the  various  potential  future  emission  scenarios.  This  work  selected  six 
representative  models  and  three  emission  scenarios  (introduced  above), 
resulting  in  18  different  forecasts.  The  International  Centre  for  Tropical 
Agriculture  (CIAT)  has  downscaled  these  climate  projections  from  the 
IPCC  GCM  model  results  (Ramirez  and  Jarvis  2010).  Their  approach, 
called  the  delta  method,  compares  the  forecasted  average  monthly 
temperature  and  rainfall  differences  for  each  GCM  grid  cell  (typically  2-3 
degrees  resolution)  with  the  historic  record  and  then  applies  that  change 
to  the  WorldClim  databases  at  30  arc-seconds  resolution.  The  result  is  a 
detailed  forecast  of  the  consequences  of  climate  change  at  less  than  a  l-km 
resolution  over  the  earth’s  landmass. 

2.7  Bioclimatic  parameters 

The  basic  WorldClim  and  CIAT  maps  provide  average  temperature  and 
rainfall  for  each  month  over  historic  and  future  time  periods.  From  this 
information,  a  set  of  secondary  datasets  are  created  to  include  20 
bioclimatic  concerns  useful  in  characterizing  the  biological  environment 
based  on  the  predicted  GCM  changes.  These  20  parameters  (Table  1) 
represent  many  of  the  concerns  that  affect  the  living  environment  in  a 
locality.  They  are  derived  directly  from  the  base  temperature  and 
precipitation  data.  Characteristics  of  the  data  are  well  documented  at  the 
CIAT  web  site.  These  are  the  data  types  used  in  the  MVA. 
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Table  1.  Bioclimatic  categories  used  for  climate  change  evaluations. 


Derived  from  maximum  and  minimum  temperature  (°C  x  10): 

BI01  =  Annual  Mean  Temperature 

BI02  =  Mean  Diurnal  Range  (Mean  of  monthly  (max  temp  -min  temp)) 

BI03  =  Isothermality  (mean  diurnal  range/temperature  annual  range) 

BI04  =  Temperature  Seasonality  (standard  deviation  *100) 

BI05  =  Max  Temperature  of  Warmest  Month 

BI06  =  Min  Temperature  of  Coldest  Month 

BI07  =  Temperature  Annual  Range  (P5-P6) 

BI08  =  Mean  Temperature  of  Wettest  Quarter 

BI09  =  Mean  Temperature  of  Driest  Quarter 

BI010  =  Mean  Temperature  of  Warmest  Quarter 

BI011  =  Mean  Temperature  of  Coldest  Quarter 

Derived  from  precipitation  (in  centimeters) : 

BI012  =  Annual  Precipitation 

BI013  =  Precipitation  of  Wettest  Month 

BI014  =  Precipitation  of  Driest  Month 

BI015  =  Precipitation  Seasonality  (Coefficient  of  Variation) 

BI016  =  Precipitation  of  Wettest  Quarter 

BI017  =  Precipitation  of  Driest  Quarter 

BI018  =  Precipitation  of  Warmest  Quarter 

BI019  =  Precipitation  of  Coldest  Quarter 

Additional  dataset: 

BI020  =  Consecutive  Months  -  the  maximum  number  of  consecutive  dry  months  of 
<100  mm  in  a  year 

2.8  Averaging  the  data 

Six  different  climate  models  and  three  different  emission  scenarios  produce 
18  sets  of  forecast  results  for  three  future  time  periods.  This  study  averaged 
the  forecasts  for  each  of  the  time  intervals.  That  is,  the  18  results  (from  the 
six  GCMs  and  three  scenarios)  were  averaged  together  for  each  of  the  maps 
in  Table  l  and  for  each  of  the  future  time  periods.  The  results  represent  the 
scientific  consensus  of  the  best  models  available.  It  also  eliminates  the  con¬ 
troversy  that  seems  to  always  emerge  when  extreme  values  from  one  model 
are  unfavorably  compared  to  those  of  another.  Finally  it  simplifies  and  clari¬ 
fies  the  task  of  comparing  the  results  of  the  MVA. 
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3  Exploration  of  the  Climate  Change  Data 

The  WorldClim  and  associated  CIAT  dataset  is  based  on  historic  rainfall 
and  temperature  readings  from  weather  stations  around  the  world.  Local 
equipment  or  settings  for  the  weather  station  can  result  in  readings  that 
are  uncoordinated  with  surrounding  weather  stations.  The  interpolation 
algorithms  used  to  populate  a  field  of  readings  between  stations  can  result 
in  slight  errors  that  can  be  revealed  as  “holes,”  circular  patterns,  or  “bulls 
eyes,”  depending  on  the  color  table  used.  The  area  near  Fort  Benning  can 
be  characterized  as  a  region  of  rolling  low  hills.  In  such  an  area,  one  would 
expect  that  the  climatic  factors  will  smoothly  vary  as  if  a  curtain  were 
draped  over  the  terrain.  Locally,  the  smoothness  of  the  curtain  would  be 
expected  to  moderate  only  because  topography  has  subtlety  been 
integrated  into  the  bioclimatic  data  through  the  WorldClim  downscaling 
process.  For  the  area  in  close  proximity  to  Fort  Benning,  the  data  show 
something  dramatically  different  than  expected. 

Examination  of  the  original  maps  revealed  an  interesting  pattern.  Almost 
centered  on  the  installation  is  a  local  data  anomaly  (Figure  l).  In  this  image, 
the  highly  enhanced  color  table  proceeds  from  red,  to  dark  red,  to  green,  to 
dark  green,  to  blue,  and  finally  to  dark  blue  as  the  precipitation  increases. 

The  climatic  pattern  does  not  just  correlate  with  topography,  which  can  be 
seen  subtlety  matching  stream  patterns.  The  Benning  bull’s  eye  is  evident 
in  the  maps  of  those  bioclimatic  factors  listed  in  Table  2. 

Table  2  lists  eight  concerns  that  clearly  exhibit  the  bull’s  eye.  Seven  are  as¬ 
sociated  with  precipitation  while  only  one  is  exclusively  related  to  temper¬ 
ature.  This  means  that  precipitation  is  the  overriding  concern  contributing 
to  the  climatic  anomaly  pattern  in  most  of  the  secondary  data  layers.  The 
image  in  Figure  3  suggests  that  Fort  Benning  is  more  intrinsically  related 
to  the  precipitation  amounts  65  km  to  the  east  of  the  installation  than  to 
areas  immediately  surrounding  the  Fort.  Anyone  studying  Fort  Benning  in 
isolation  or  the  immediately  adjacent  area  would  not  detect  this  feature 
since  its  presence  does  not  emerge  until  the  study  area  is  drawn  out  to 
roughly  15  km  beyond  the  installation’s  eastern  boundaries. 
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Figure  1.  Biol2  annual  precipitation  exhibits  the  Benning  “climatic  hole".  The  darker  blue 
represents  locally  lower  precipitation.  The  image  is  in  units  of  standard  deviations  from  the 
dataset  mean  (explained  later).  The  color  table  has  been  stretched  greatly  to  enhance  the 
"hole."  Xs  identify  the  locations  of  permanent  weather  stations. 


Table  2.  Appearance  of  the  Benning  bull’s  eye  in  the  bioclimatic  layers. 


BioNumber 

Meaning 

Associated  with 
Precipitation 

Associated  with 
Temperature 

Bio8 

Mean  Temperature  of  Wettest  Quarter 

Yes 

Yes 

BiolO 

Mean  Temperature  of  Warmest  Quarter 

No 

Yes 

Biol2 

Annual  Precipitation 

Yes 

No 

Biol3 

Precipitation  of  Wettest  Month 

Yes 

No 

Biol4 

Precipitation  of  Driest  Month 

Yes 

No 

Biol6 

Precipitation  of  Wettest  Quarter 

Yes 

No 

Biol7 

Precipitation  of  Driest  Quarter 

Yes 

No 

Biol9 

Precipitation  of  Coldest  Quarter 

Yes 

Yes 

The  hole  can  be  thought  of  as  an  anomaly  in  the  smoothly  draped  curtain  of 
data.  In  fact,  it  is  a  thumbprint  depression,  which  in  essence  extends  the 
precipitation  character  that  exists  to  the  east  pulling  it  toward  the  location 
of  Fort  Benning.  It  is  significant  in  that  it  sits  over  the  target  Army 
installation  like  a  “bulls  eye”  mark.  It  suggests  that  an  errant  set  of  readings 
may  originate  from  the  southern  quarter  of  the  installation.  In  Figure  l, 
each  weather  station  is  labeled  by  name  plotted  with  an  “x.”  The  weather 
stations  are  the  source  for  the  data  that  the  WorldClim-based  climate 
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downscaling  process  uses.  In  fact,  the  Fort  Benning  weather  station  does 
reside  in  the  hole,  though  not  at  its  center.  It  must  be  suspected  that 
readings  from  this  station  might  be  the  source  of  the  data  hole. 

Figure  2  displays  the  CIAT  Bioi2  annual  precipitation  map.  Sixteen  or 
more  bull’s  eyes  can  be  seen.  Evidently,  when  the  data  was  being 
generated,  the  smoothing  algorithm  spawned  local  holes  and  hills.  Since  it 
seems  to  be  a  common  occurrence,  it  is  apparently  one  of  the  drawbacks  of 
using  downscaled  data. 

The  point  of  this  discussion  is  that  this  bull’s  eye  artifact  has  multiplied 
itself  in  some  of  the  derived  data  layers.  It  is  therefore  predicted  that  the 
Benning  Hole  will  appear  in  the  results  of  the  multivariate  classifications. 
Also,  it  can  be  expected  that  patches  50  km  to  the  east  of  the  installation  will 
be  more  similar  to  the  installation  than  will  areas  either  on  the  installation 
or  more  immediately  adjacent  to  it.  In  fact,  many  of  the  resulting  category 
identifications  do  have  a  very  strong  east-west  bias.  Specifically,  it  appears 
in  both  the  Unsupervised  and  Supervised  Classification  results. 

Figure  2.  Multiple  sample  "bull’s  eyes”  visible  in  the  Annual  Precipitation  data.  Units  are  in 
standard  deviations  from  the  image  mean.  Deep  blue  is  the  lowest  precipitation,  white  is  near 
the  mean,  and  deeper  reds  are  the  highest  values.  Circles  show  obvious  holes.  Different  color 
tables  enhance  different  anomalies.  Not  all  holes  coordinate  with  weather  stations. 
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4  Multivariate  Analyses 

The  goal  of  this  work  is  to  identify  current  locations  that  have  the 
combination  of  climate  characteristics  that  are  anticipated  for  the 
installation,  Fort  Benning,  in  the  future.  One  approach  is  to  simply  identify 
areas  that  have  exactly  the  same  set  of  characteristics.  Unfortunately,  there 
is  no  guarantee  that  a  current  location  will  have  the  identical  combination 
of  values  anticipated  in  the  future.  By  setting  acceptable  ranges  of  values 
for  each  of  the  variables,  it  is  possible  to  begin  to  identify  similar  areas. 

The  challenge  that  immediately  arises  is  how  to  establish  ranges  that  are 
equally  valid.  The  same  challenge  exists  in  the  processing  of  satellite  (or 
high  altitude)  imagery  to  identify  locations  of  similar  land  cover.  Such 
imagery  is  typically  associated  with  grey-scale  images  representing  several 
to  many  different  wavelengths  that  are  statistically  processed  using  MVA, 
which  has  been  adopted  to  address  the  challenge  of  locating  similar  areas. 

4.1  Background  for  multivariate  analyses 

An  MVA  is  based  on  the  statistical  principle  that  involves  observation  and 
analysis  of  more  than  one  statistical  outcome  variable  at  a  time.  In  the  area 
of  land  planning,  MVAs  have  been  traditionally  applied  to  satellite  images. 
In  this  area,  the  variables  are  normally  seen  as  the  three  grey-scale  images 
taken  by  a  spectral  sensor  in  the  red,  green,  blue,  and  often  infrared  parts 
of  the  spectrum.  However,  multispectral  satellites  often  sense  five  to  about 
15  distinct  spectral  bands.  Newer  satellite  imaging  systems  generate  many 
images  each  taken  in  a  narrow  band  width  in  the  radiation  spectrum. 

These  newer  satellite  imaging  systems  are  called  hyper-spectral  systems. 
Traditionally  values  from  these  bands  have  been  submitted  to  a  mathe¬ 
matical  manipulation  routine  (an  MVA  program)  to  generate  patches  that 
represent  distinct  concerns,  normally  those  concerns  being  distinct  land 
cover  types.  Mathematicians  have  developed  several  distinct  algorithms  to 
analyze  the  input  bands.  These  multivariate  programs  are  now  readily 
available  for  application  in  software  packages.  With  the  increase  in  the 
number  of  spectral  bands  available,  programs  have  been  written  to  handle 
a  large  number  of  bands  at  the  same  time.  However  these  programs  do  not 
distinguish  whether  an  input  band  is  a  visual  “photographic”  style  satellite 
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image  or  set  of  data,  as  long  as  there  is  an  X  and  a  Y  component  on  a  com¬ 
patible  Cartesian  system  and  a  value  (also  called  a  Z)  at  a  particular  point. 
The  obvious  example  of  a  non-image  dataset  is  that  of  topographic  eleva¬ 
tion  where  the  z-values  really  do  physically  represent  the  third  dimension 
of  height.  Thus,  an  elevation  map  can  be  considered  another  image, 
though  it  is  by  no  means  a  spectral  image. 

It  is  easy  to  extend  the  concept  to  allow  the  submittal  of  red,  green,  blue 
imagery  bands  along  with  an  elevation  dataset  to  a  multivariate 
classification  program.  The  result  of  this  analysis  might  therefore  be  land 
uses  distinguished  by  elevation,  as  well  as  spectral  reflection.  Multivariate 
images  arise  from  a  surprising  variety  of  sources:  temperature, 
gravitational  field,  impedance,  magnetic  field,  electrical  field,  mass, 
wavelength,  ultrasound  wavelength,  polarization,  electron  energy,  etc. 
Common  subdivisions  of  the  scientific  imaging  fields  includes  satellite 
imaging,  medical  (clinical)  imaging  and  in  microscopy.  As  long  as  an 
image  has  two  pixel  indices  (e.g.,  an  X  and  Y  location  in  the  image  plane) 
and  a  variable  quantity  to  make  a  three-dimensional  array,  it  can  be  used 
as  one  of  the  components  in  an  MVA  analysis. 

For  the  purposes  of  this  work,  it  is  desirable  to  substitute  climatic  indices 
for  spectral  bands  for  submission  to  MVA  routines  to  determine  how 
climate  change  will  alter  the  environment  around  military  installations 
and  to  predict  what  the  new  climatic  conditions  will  be  in  the  future.  Fort 
Benning  GA  is  the  focus  of  the  actual  analyses  here,  though  they  can  be 
applied  to  any  other  installation. 

The  application  of  MVA  techniques  using  climate  data  has  been 
demonstrated  (Hargrove  2004,  Westervelt  2011).  In  this  work,  four 
techniques  were  successfully  applied:  (1)  Unsupervised  Classification, 

(2)  Supervised  Classification,  (3)  Principal  Components,  and  (4)  the 
Primarily  Analogous  Multivariate  approach.  Canonical  Components 
Analysis  (CCA)  and  Sequential  Maximum  A  Posteriori  (SMAP)  estimation 
were  attempted,  but  no  usable  signature  file  could  be  generated.  Others 
like  Tassel  Cap  are  specific  to  particular  satellite  sensors  and  were  not 
appropriate  for  this  application  (ERDAS  2010). 
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4.2  Transforming  the  bioclimatic  indices  for  MVA  analysis 

The  bioclimatic  maps  needed  to  be  prepared  for  submission  to  the 
multivariate  classifiers.  The  multivariate  classifiers  usually  expect  the  data 
to  have  the  following  characteristics: 

•  They  must  range  between  o  to  255  (inclusive). 

•  The  histogram  distribution  should  have  a  semi-normal  distribution 
(i.e.,  tails  that  contain  few  values  with  values  peaking  at  one  value). 

•  The  histogram  distribution  should  be  continuous  . 

Part  of  the  problem  with  the  bioclimatic  data  layers  is  that  ranges  can  be 
negative  (e.g.,  temperatures  below  o  °C);  an  index  range  can  vary  from  -1 
to  +1,  and  maximum  values  can  range  into  the  thousands.  To  make  the 
different  data  layers  compatible,  it  was  necessary  to  normalize  the  data 
without  losing  the  information  value  it  contained. 

One  might  think  that  an  obvious  and  straightforward  approach  would 
simply  involve  scaling  data  values  from  their  actual  minima  and  maxima  to 
a  scale  from  o  to  255.  Unfortunately  this  is  not  appropriate  because  it 
assumes  that  the  ranges  and  increments  are  equally  important.  For 
example,  if  the  range  of  the  raw  bioclimatic  parameter  of  X  for  2025  was  25 
to  too  and  the  range  for  the  comparable  2055  data  was  25  to  80,  when  these 
layers  were  rescaled  between  o  to  255,  both  would  range  from  o  to  255, 
implying  the  maximum  values  were  the  same.  However  it  is  clear  that  the 
data  for  2055  show  an  important  change  -  a  decrease  in  the  maximum.  This 
important  change  is  lost  in  rescaling  from  o  to  255.  A  means  of  normalizing 
all  the  data  layers  while  preserving  the  character  of  the  raw  data  is  needed. 

The  technique  adopted  was  to  generate  layers  that  showed  how  far  each 
value  was  from  the  layer  mean  as  measured  in  units  of  standard  deviations. 
If  the  full  data  set  is  used  to  generate  the  resultant  values,  these  values  are 
called  the  “z  values.”  If  a  sampling  of  all  the  data  is  used  to  generate  the 
resultant  values,  the  values  are  call  “students-t  distribution.”  Both  are 
common  statistical  techniques  (Sprinthall  2003).  The  technique  used  is  a 
cross  between  the  two.  More  importantly,  the  results  are  all  in  the  same 
units,  standard  deviations. 
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The  first  step  is  to  calculate  the  mean  and  standard  deviation  across  the 
population  for  each  layer.  For  each  bioclimatic  parameter,  the  population 
included  the  baseline  data  AND  all  the  values  for  2025,  2055,  and  2085. 
The  next  step  is  to  return  to  each  of  the  original  time  horizon  maps  for 
each  bioclimatic  parameter  and  calculate  the  z-value  for  each  location, 
done  here  with  the  (simple)  Equation  1: 

BioclimN(date)z  =  (bioclirriN(date)value  -  bioclimN(all)avg)  /bioclimN(all)sd  Eq  1 


where: 


BioclimN 

Date 

z 

value 

(all)avg 

(all)sd 


the  Nth  Bioclimate  concern  of  the  19 

1990,  2025,  2055  or  2085 

the  z-value  for  that  grid  cell 

the  input  value  of  that  grid  cell 

the  average  for  that  Bioclim  layer  that  includes  data 

from  all  four  time  horizons 

the  Standard  Deviation  for  that  Bioclim  layer  that 
includes  data  from  all  four  time  horizons. 


An  increment  of  127  was  added  to  make  this  value  better  fit  in  the  required 
0-255  range,  to  ensure  no  negative  values  existed  and  to  make  all  the  means 
compatible.  The  final  ranges  almost  always  fell  into  the  range  of  124  to  130. 
Thus  the  final  maps  almost  always  fulfilled  all  three  requirements.*  Tests 
were  conducted  to  see  if  stretching  the  124  to  130  range  to  1  to  255  would 
make  any  difference  in  the  final  results  from  the  Unsupervised 
Classification  program.  The  test  results  were  the  same  with  or  without  the 
stretch  so  the  unstretched  data  were  used  for  the  rest  of  the  analyses 

Where  the  range  of  raw  values  was  very  small,  the  resultant  map  might 
have  a  conveniently  normal  distribution,  but  it  was  often  discontinuous, 
thus  failing  the  third  requirement.  Examples  of  bioclimatic  maps  that 
failed  the  third  requirement  are: 

•  Bio3  Isothermality  has  only  18  levels 

•  Consecutive  Dry  Months  has  only  12  levels. 


*  It  is  possible  but  highly  unlikely  a  value  below  zero  or  above  255  may  exist  in  the  resultant  data.  To 
prevent  these  unlikely  events,  all  values  above  255  were  assigned  to  255;  all  values  below  zero  were 
assigned  to  zero.  These  extremes  however  never  occurred  in  the  output  data 
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Thus  the  use  of  Bio3  and  Cons_mths  were  limited  because  they  met  the 
expectations  of  a  multivariate  input  very  poorly.  Tests  indicated  that  the 
classification  programs  would  still  work  in  spite  of  this  characteristic.  This 
was  fortunate  because  it  was  suspected  that  consecutive  dry  months  might 
become  an  important  ecosystem  distinguishing  factor. 

4.3  Applying  Unsupervised  Classification 

Unsupervised  classification  is  a  computer-automated  MVA  method  of 
pattern  recognition  in  which  parameters  specified  by  the  user  are  used  to 
group  data  into  statistically  separated  sets  of  locations.  The  result  is  a  set 
of  statistically  distinct  data  “signatures”  that  can  be  used  to  categorize  each 
location,  which  often  results  in  the  creation  of  a  map  of  statistically 
distinct  patches.  This  is  because  typically  neighboring  locations  share 
similar  signatures.  Carrying  out  an  Unsupervised  Classification  consists  of 
two  steps:  clustering  all  the  data  into  statistically  separate  group 
identifications  and  then  assigning  each  grid  cell  to  the  signature  that  is 
most  similar  to  their  spectral  characteristics. 

The  clustering  algorithm  (GRASS  Development  Team  2012)  reads  through 
the  raster  data  and  builds  sets  of  locations  that  share  a  similar  spectral 
reflectance  (see  Figure  3).  The  pixel  clusters  are  imagery  categories  that 
can  be  related  to  land  cover  types  on  the  ground.  The  spectral  distributions 
of  the  clusters  (which  will  be  the  land  cover  spectral  signatures)  are 
influenced  by  parameters  set  by  the  user.  Parameters  include  number  of 
resultant  classes,  sample  interval,  iterations  limit,  convergence,  cluster 
separation,  and  minimum  number  of  pixels  for  a  valid  class.  In  effect, 
when  parameters  were  set  they  always  tended  to  decrease  the  number  of 
resultant  classes.  This  work  strived  to  get  as  many  classes  as  possible. 
Therefore  the  parameters  were  almost  always  the  default  values. 

The  clustering  starts  by  generating  spectral  signatures  for  a  number  of 
clusters  and  by  “attempting”  to  end  up  with  this  number  of  clusters  during 
the  clustering  process.  The  cluster  spectral  signatures  that  result  are 
composed  of  cluster  means  and  covariance  matrices.  The  signature  file  used 
as  input  for  the  second  step  is  the  file  generated  by  the  clustering  routine. 
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Figure  3.  The  clustering  algorithm  builds  pixel  clusters  (right)  based  on  the  spectral 

reflectance’s  of  the  pixels. 


In  the  second  step,  these  cluster  means  and  covariance  matrices  are  used 
to  classify  all  the  pixels  by  separating  them  into  groups  on  the  basis  of  the 
maximum  likelihood  that  the  individual  pixel  belongs  in  a  specific  group. 
In  this  case,  the  maximum  likelihood  is  equal  to  the  minimum  statistical 
distance  from  the  mean  of  any  of  the  groups  possible.  The  clusters  or 
spectral  classes  that  result  can  be  related  to  land  cover  types  on  the  ground 
(Figure  3,  right).  The  user  has  to  specify  the  name  of  a  group,  the  name  of 
subgroup,  and  the  name  of  a  file  that  contains  the  signatures.  The 
subgroup  contains  only  the  imagery  band  files  that  the  user  wishes  to 
classify.  The  discussion  below  compares  the  effect  that  using  a  different 
number  of  bands  has  on  the  output  resultant  map.  This  research 
substituted  bioclimatic  bands  for  the  traditional  imagery  bands. 

As  alluded  to  previously,  in  the  analysis  of  the  bioclimatic  data  layers 
transformed  into  z-values,  signatures  for  clusters  were  developed  by  sim¬ 
ultaneously  analyzing  1990,  2025,  2055,  and  2085  datasets.  These  clusters 
were  then  used  to  categorize  each  set  independently  to  create  a  single  cat¬ 
egorized  map  for  each  of  the  four  time  frames.  This  meant  that  the  Catego¬ 
ry  1  in  each  map  is  associated  with  the  same  signature,  allowing  identifica¬ 
tion  of  areas  that  exist  today  that  share  characteristics  of  the  future  Fort 
Benning.  Running  the  Unsupervised  Classification  on  the  data  gave  com¬ 
parable  results  for  the  1990  and  2025  time  frames.  Analyses  were  run  for 
2055  and  2085,  but  resulted  in  no  patches  that  were  the  same  as  those  at 
Fort  Benning  in  1990  so  no  comparisons  for  these  dates  can  be  carried  out. 
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4.3.1  Effects  of  variable  layer  inputs  on  unsupervised  classifications  and 
class  migration 

This  section  integrates  the  discussion  of  the  results  into  a  description  of 
the  effect  of  varying  the  number  of  input  bands  or  layers  in  what  is 
basically  a  sensitivity  analysis.  The  number  of  bands  is  being  varied  to 
determine  if  adding  more  bands  produces  a  better  output. 

Varying  the  number  of  input  bands  results  in  a  problem  in  identifying  the 
categories  that  are  the  same  between  submissions.  To  simplify  this  and 
make  the  discussion  relevant,  this  work  focused  on  what  occurred  to  those 
classified  patches  that  best  characterized  the  immediate  Fort  Benning  area. 
It  was  assumed  that  a  patch  located  in  the  same  place  at  the  same  time  in 
each  of  the  three  trials  is  equivalent  between  analyses  with  a  different 
number  of  input  layers.  The  research  question  of  interest  then  becomes: 

From  where  in  1990  did  the  patch  that  exists  at  Fort  Benning  for  the 
period  2025  come? 

Using  a  sample  of  30  patches  in  each  time  horizon  (1990  and  2025)  and 
each  multilayer  submission  (4,  6,  and  8),  the  characteristics  of  patch 
movement  like  those  in  Figure  4  were  measured.  Table  3  lists  the  averaged 
values  and  the  results. 

Movement  in  the  X  (east-west)  direction  is  small.  The  movement  in  the  Y 
(north-south)  direction  slowly  decreases  as  more  layers  are  added  and 
always  varies  more  than  that  in  the  X  direction. 

For  four  layers,  two  categories  dominate  the  region  around  Fort  Benning 
in  2025.  Category  19  comes  from  a  location  in  1990  (in  Figure  4  Category 
19  is  light  blue)  that  is  about  130  km  to  the  southeast  of  the  installation 
while  Category  75  on  Fort  Benning  comes  from  an  original  patch  (the 
lighter  green)  that  was  120  km  to  the  south-southeast  of  the  installation. 


Table  3.  Cluster  migration  amounts  1990-2025. 


For  4  Layers 

Bios  1,  9, 12, 19 

For  6  Layers 

Bios  1,  9, 11, 12, 14, 19 

For  8  Layers 

Bios  1,  2, 9, 11, 12, 13, 14, 19 

Change  in 

X 

y 

X 

y 

X 

y 

Average* 

5.8 

12.1 

-3.0 

7.9 

0.1 

5.1 

SD 

10.7 

9.1 

16.0 

7.2 

7.7 

5.0 

‘Values  are  unitless. 
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For  six  layers,  Category  35  dominates  southwestern  Fort  Benning  and 
Category  30  dominates  northeastern  Fort  Benning  in  2025.  Category  30 
comes  from  a  location  in  1990  that  is  about  55  km  to  the  south  of  the 
installation.  (In  Figure  5,  Category  30  in  1990  is  light  green.)  Category  35  on 
Fort  Benning  comes  from  an  original  patch  (the  lighter  blue)  that  was  140 
km  to  the  southeast  of  the  installation.  This  is  the  effect  of  the  “Benning 
hole”  as  predicted. 

For  eight  layers  (Figure  6),  Category  18  dominates  Fort  Benning  in  2025. 
Category  18  comes  from  a  location  in  1990  that  is  nearly  too  km  to  the 
south-south  east  of  the  installation. 

In  all  three  cases,  even  though  a  different  number  of  input  layers  were 
used,  the  source  of  the  Fort  Benning  bioclimatic  characterization  for  2025 
originates  in  1990  from  about  too  km  to  the  southeast.  This  implies  that, 
in  2025,  the  Fort  Benning  natural  environment  will  become  similar  to  the 
1990  environment  in  Sasser,  GA. 

4.3.2  Ecosystem  migration  to  Fort  Benning 

Overlaying  a  map  of  ecosystems  onto  this  region  shows  how  the 
ecosystems  will  migrate  in  the  Fort  Benning  region  (Figure  7).  The 
installation  currently  lies  at  the  junction  of  two  major  eco-units.  Part  of 
Fort  Benning  (the  yellow  southwest  third)  is  in  the  Coastal  Plains  -  Middle 
Section.  The  rest  of  it  is  in  the  Gulf  Coastal  Plains  and  Flatwoods  Section. 
Figure  7  shows  that,  from  the  Forest  Service  Ecosystem  delineation,  the 
location  too  km  to  the  south-southeast  that  will  migrate  to  the  Fort 
Benning  area  by  around  2025  will  be  a  combination  of  the  Wiregrass 
Plains  and  the  Plains  Subsections.  Although  most  of  Fort  Benning  is 
currently  at  the  northern  edge  of  the  Gulf  Coastal  Plains  in  the  Sandhills 
Subsection,  this  does  not  mean  that  its  ecosystem  will  be  changed  by  2025. 
There  will  be  a  lag  in  the  occupying  vegetation. 
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Figure  4.  Four  layer  clustering  results  for  Fort  Benningfor  1990-2025. 


Figure  5.  Six  layer  clustering  results  for  Fort  Benningfor  1990-2025. 
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Figure  6.  Eight  layer  clustering  results  for  Fort  Benning  for  1990-2025 


Figure  7.  Migration  of  ecosystems  into  the  Fort  Benning  region  based  on  the  Unsupervised 

Classification  results  for  2025. 


ERDC/CERL  TR-15-2 


21 


However,  it  does  mean  that  those  climatic  parameters  that  helped  to  form 
the  Sandhills  will  have  migrated  away  to  be  replaced  by  new  parameters 
that  are  characteristic  of  The  Plains  and  Wiregrass  Plains  areas.  That  is, 
Fort  Benning’s  ecosystem  will  attempt  to  adapt  to  climatic  forces  formerly 
associated  with  the  Wiregrass  Plains  and  the  Plains.  Vegetation  and  maybe 
fauna  will  be  influenced.  Of  course  the  sandy  soils  of  the  Sandhills  will  not 
migrate.  This  observation  suggests  that  a  new  nomenclature  for 
ecosystems  must  be  developed  if  researchers  are  actually  going  to  follow 
their  changes  in  over  time. 

Table  4  lists  some  of  the  differences  between  the  current  and  future  Fort 
Benning  ecosystems.  The  data  in  Table  4  make  it  clear  that  not  only  will 
the  ecosystems  change,  but  some  of  the  important  ecosystem 
characteristics  will  alter.  A  major  change  is  that  expected  snowfall  will 
disappear.  Snow  reflects  a  great  deal  of  light  back  into  space.  Without  the 
snow,  that  radiation  will  be  absorbed.  This  will  contribute  to  winter  time 
heating  significantly.  Along  with  this,  the  mean  monthly  minimum 
temperature  in  January  will  rise  from  below  freezing  (320  F)  to  above 
freezing  by  several  degrees.  This  is  important  because  many  plant  and 
animal  species  cannot  survive  freezing  temperatures.  One  can  expect  many 
more  types  of  invasive  species  to  be  able  to  establish  a  foothold  in  the  Fort 
Benning  area  than  would  be  expected  due  to  a  simple  rise  in  temperature 
of  a  few  degrees.  The  number  of  frost-free  days  will  increase  significantly 
so  the  growing  season  will  be  about  a  month  longer.  The  mean 
temperature  will  increase  by  over  1  degree  Fahrenheit  and  the  summer 
humidity  will  also  increase  by  2-3%.  Fort  Benning  will  be  warmer  in  the 
winter  and  even  more  hot  and  humid  in  the  summer. 

It  is  useful  to  point  out  here  that  the  preceding  description  is  not  based  on 
the  bioclimatic  data  itself.  Instead  the  data  determined  the  multivariate 
categories  that  were  then  used  to  find  the  source  location  (and  therefore 
ecosystem)  that  would  begin  to  establish  itself  at  Fort  Benning  in  2025. 

The  description  of  the  ecosystem  at  that  location  in  1990  was  then  used  to 
characterize  the  Fort  Benning  area  in  2025.  This  is  a  valid  but  less  direct 
approach  to  arrive  at  the  same  information  that  could  be  extracted  from 
the  original  bioclimatic  data. 
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One  significant  observation  is  that  Fort  Benning  resides  on  the  northern 
edge  of  a  major  ecosystem  division,  the  Gulf  Coastal  Plains  and  Flatwoods 
Section  (versus  the  Coast  Plains  Middle  Section).  Species  near  the  north¬ 
ern  boundaries  tend  to  be  “hardier”  then  their  relatives  in  the  middle  of 
their  range  or  near  the  southern  edge.  As  Fort  Benning’s  bioclimatic  char¬ 
acteristics  change  to  what  is  now  a  more  southerly  eco-type,  those  species 
present  will  be  less  stressed  so  their  hardiness  can  be  expected  to  decrease. 

4.4  The  Supervised  Classification  approach 

In  the  Supervised  Classification  approach,  the  signatures  associated  with 
potential  classes  are  based  on  user  identification  of  the  locations  of  certain 
areas.  When  used  in  association  with  image  processing,  the  user  defines 
areas  by  drawing  boundaries  around  the  limits  of  important  features.  For 
example,  in  a  satellite  image,  the  user  may  want  to  distinguish  between 
soybean  and  cornfields  so  he/she  defines  a  few  areas  that  contain  each 
type.  The  classifying  routine  finds  the  spectral  signatures  of  each  type  and 
uses  this  characterization  to  find  other  similar  areas.  In  this  case,  it  was 
desired  to  define  an  area  (Fort  Benning)  and  find  similar  areas  over  time 
rather  than  within  the  same  image. 

To  carry  out  this  procedure,  it  was  necessary  to  first  define  sample  locations 
representing  areas  of  interest  in  the  form  of  a  training  map.  In  this  case,  the 
training  map  was  the  topographic  elevation  map  reclassified  into  57 
elevation  classes.  Many  of  the  climatic  concerns  mimic  and/or  are 
influenced  by  topography  so  using  topography  was  a  means  to  generate 
marginally  related  but  independent  analysis  for  the  training  map  input.  A 
signature  was  then  calculated  for  each  of  these  areas,  consisting  of  value 
mean  and  standard  deviations  for  each  of  seven  bioclimatic  maps  (Table  5). 


Table  5.  Bioclimatic  characteristics  used  in  the  Supervised  Classification  analysis. 


Layer 

Characteristics 

Bioclimatic  Layer  1 

Annual  Mean  Temperature 

Bioclimatic  Layer  2 

Mean  Diurnal  Range  (Mean  of  monthly  (max  temp  -min  temp)) 

Bioclimatic  Layer  9 

Mean  Temperature  of  Driest  Quarter 

Bioclimatic  Layer  11 

Mean  Temperature  of  Coldest  Quarter 

Bioclimatic  Layer  12 

Annual  Precipitation 

Bioclimatic  Layer  13 

Precipitation  of  Wettest  Month 

Bioclimatic  Layer  19 

Precipitation  of  Coldest  Quarter 
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Notice  that  these  are  not  the  same  layers  as  used  previously.  The  original 
intent  was  to  always  use  the  same  set  of  bioclimatic  layers  as  the  inputs 
between  the  analyses.  It  was  found  that  success  in  getting  any  output 
varied  depending  on  which  layers  were  used.  (Also  the  number  of  output 
categories  varied  greatly  depending  on  the  combinations  of  inputs  used.) 
This  research  presents  the  best  outputs  obtained.  Usually,  a  large  number 
of  input  combinations  were  tested  in  the  process.  Even  so,  as  is  evident 
from  the  discussions,  it  was  often  possible  to  follow  a  specific  category 
across  only  a  couple  time  horizons  rather  than  all  four. 

At  this  point,  it  is  desirable  to  identify  the  current  location  of  the  areas  that 
most  closely  match  the  conditions  anticipated  at  Fort  Benning  in  2085. 

The  established  signatures  were  applied  to  the  set  of  maps  associated  with 
2085  using  a  maximum  likelihood  classifier  and  it  was  discovered  that 
Fort  Benning  is  associated  with  a  single  signature  (Category  16  shown  in 
Figure  8a).  Applying  the  signature  set  to  the  suite  of  2055  (Figure  8b), 
2025  (Figure  8c),  and  current  maps  (Figure  8d)  allows  identification  of  the 
areas  that  are  assigned  to  Category  16.  The  last  lower  right  quadrant  of  the 
image  identifies  where  one  might  go  today  to  experience  the  bioclimatic 
conditions  anticipated  in  2085. 

Figure  8  shows  that  there  is  a  march  of  similar  areas  (in  blue)  moving  to 
the  north  to  Fort  Benning.  In  fact,  it  is  easier  to  understand  these  results  in 
reverse  order  -  in  the  last  image  (1990)  the  2085  climate  exists  just  north 
of  the  Florida  Panhandle  state  line.  As  time  passes,  the  areas  similar  to  the 
2085  characterization  slowly  march  roughly  north  and  west  to  finally  be  at 
Fort  Benning’s  location  by  2085.  The  analysis  results  for  2025  are  similar 
to  those  previously  found  using  the  Unsupervised  Classification  approach 
described  in  Section  4.3.1.  (The  unsupervised  analysis  produced  no  results 
for  2055  or  2085  that  could  be  used  for  comparison  with  these.) 

4.5  The  principal  components  analysis 

Principal  Components  Analysis  (PCA)  is  a  method  of  data  compression 
that  allows  less  redundant  data  to  be  expressed  in  lower  band  numbers. 

For  example,  if  a  user  inputs  four  bands  (normally  spectral  imagery 
bands),  the  classification  routine  will  output  four  new  bands. 
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Figure  8.  Supervised  classification  Category  16  that  represents  Fort  Benning  in  the  2085 
period  (upper  left),  the  2055  period  (upper  right),  the  2025  period  (lower  left),  and  the  1990 

period  (lower  right). 


These  new  bands  will  have  their  axes  mathematically  twisted  so  a  scatter 
plot  of  the  bands  of  data  represent  the  widest  variance  among  the  data  in 
the  first  band  and  each  successive  axis  (band)  represents  the  next  smaller 
amount  of  variance  that  is  not  already  represented.  Principal  components 
axes  are  orthogonal  (perpendicular)  to  one  another.  Principal  component 
bands  are  uncorrelated  and  nonredundant  since  each  principal  component 
describes  different  variance  within  the  original  data. 

The  drawback  in  using  the  traditional  Unsupervised  Classification 
technique  is  that,  if  the  underlying  variance  matrices  were  poorly  formed, 
then  the  maximum  likelihood  classifier  would  not  run  or  it  would  generate 
only  a  few  classes.  Furthermore,  as  has  been  mentioned,  some  of  the 
bioclimatic  datasets  seem  to  be  highly  related  so  that  those  that  are 
“repeated”  are  given  extra  weight  in  the  final  analysis.  A  PCA  realigns  the 
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data  axes  so  that  most  of  the  variance  in  the  data  is  along  the  first  axis,  the 
next  axis  is  oriented  orthogonally  along  the  direction  of  the  next  highest 
variance  direction,  etc.  for  all  further  axes  generated.  This  means  that 
most  of  the  data  information  is  located  in  the  first  few  layers  generated 
and  that  the  variance  matrices  are  well  formed  for  use  by  the  maximum 
likelihood  classifier.  PCA  groups  are  created  by  forming  composite  axes 
that  maximize  the  overall  distance  between  the  data.  In  other  words,  PCA 
determines  the  net  effect  of  each  variable  on  the  total  variance  of  the 
dataset,  and  then  extracts  the  maximum  variance  possible  from  the  data. 
Normality  in  the  distribution  of  variables  is  not  strictly  required,  but  it 
does  enhance  the  analysis.  It  was  decided  to  apply  the  PCA  to  the  issue  of 
climate  change  at  Fort  Benning. 

As  a  quick  check  of  these  assumptions,  a  small  test  study  area  slightly 
larger  than  Fort  Benning  was  selected.  The  basic  procedure  was  applied 
using  four  precipitation  bioclimatic  layers  and  four  temperature 
bioclimatic  layers  of  the  1990  bioclimatic  data.  The  resulting  PCA  output 
had  eight  layers.  Figure  9  shows  the  first  four  components  of  the  eight 
generated.  Figure  9a  shows  the  PCA  output  axis  1  (67%  of  the  variance), 
which  is  an  integration  of  topography  along  with  precipitation.  The  two 
brightest  areas  (upper  center  left  and  right  edge)  coordinate  with 
precipitation  low  districts.  These  can  have  influence  in  the  classification 
results.  Figure  9b  shows  Axis  2  (22%  of  the  variance),  in  which  a 
progression  can  be  seen  from  upper  left  to  lower  middle  with  topography 
barely  indicated.  Figure  9c  shows  Axis  3  (6%  of  the  variance),  which  is  a 
very  similar  progression  without  the  topography.  Figure  9d  shows  Axis  4 
(3%  of  the  variance),  which  exhibits  a  hyperbolic  distribution  not  clearly 
attributable  to  any  characteristic.  These  illustrate  both  the  PCA  method  of 
handling  data  and  some  of  the  attributes  inherent  within  the  data  itself. 

All  eight  PCA  output  layers  were  run  through  the  clustering  and  the 
maximum  likelihood  classifier.  The  all  8-layer  analysis  generated  147 
classes  (100%  of  data  value  used).  When  the  analysis  was  run  with  only  the 
first  four  layers  (which  contained  98%  of  the  variance)  all  of  the  PCA 
layers  were  found  to  be  acceptable  to  the  maximum  likelihood  classifier. 
The  4-layer  analysis  generated  153  classes.  In  all  cases  the  distribution  of 
patches  is  similar  except  that  the  4-layer  analysis  was  better  defined  while 
the  8-layer  analysis  was  more  “speckled.” 
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Figure  9.  First  four  components  of  the  eight  generated  by  the  PCA  for  1990.  Axis  1  is  upper 
left,  Axis  2  is  upper  right,  Axix  3  is  lower  left,  Axis  4  is  lower  right 


Two  more  analyses  were  done.  The  3-layer  classification  (containing  95% 
of  the  variance)  generated  154  classes.  The  2-layer  classification 
(containing  90%  of  the  variance)  generated  150  classes.  The  2-layer 
analysis  was  strikingly  similar  to  the  yearly  precipitation  data  hole  shown 
in  Figure  1  (p  9).  As  more  data  are  added,  the  domination  of  precipitation 
becomes  less  noticeable.  Interesting,  the  borders  of  the  classes  generated 
in  the  outputs  only  mildly  suggest  any  topographic  influence  although  PCA 
output  Layer  1  above  (Figure  9b)  seems  to  be  predominately  topography. 
One  interpretation  of  this  is  that  phenomenon  is  that  adding  further 
increments  of  less  important  information  to  the  analysis  only  tends  to 
confuse  the  result.  The  results  of  this  test  suggested  that  carrying  out  the 
full  analysis  would  result  in  an  improved  output. 

Running  the  PCA  on  the  entire  study  area  for  all  19  bioclimatic  layers  for 
all  four  time  periods  resulted  in  a  transformation  in  which  the  first  two 
layers  accounted  for  close  to  100%  of  the  variance  in  the  original  data. 
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Thus  for  the  next  step,  the  clustering  routine,  it  was  only  necessary  to  use 
these  two  bands  to  account  for  nearly  all  the  information  value. 

After  applying  the  Unsupervised  Classification  procedure  (the  clustering 
and  maximum  likelihood  analyses)  to  the  PCA  bands,  91  categories  were 
defined  in  the  output  results.  Figure  10  shows  those  classes  that  occurred 
within  Fort  Benning  in  the  1990  time  horizon.  As  can  be  seen,  Fort 
Benning  is  located  at  the  western  end  of  a  band  of  similar  classes  that 
stretch  to  the  east  and  then  north  parallel  to  the  Atlantic  coast.  In  fact,  this 
area  coordinates  closely  with  the  Forest  Services’  ecosystem  delineation  of 
the  northern  edge  of  the  Southern  Atlantic  Coastal  Plains  and  Flatwoods 
Section  within  the  Province  type  of  Outer  Coastal  Plain  Mixed  Forest. 

Figure  11  shows  those  patches  that  by  2025  will  move  to  Fort  Benning. 
Again  the  darker  areas  show  greater  occurrence  on  the  installation.  It  is 
notable  that  the  darkest  areas  (those  most  common  on  Fort  Benning  in 
2025)  are  tending  to  move  further  from  the  east  by  about  too  km.  Just  a 
small  amount  of  these  classes  still  exist  on  the  installation. 

Figure  10.  PCA  resultant  classes  that  occurred  within  Fort  Benning  in  the  1990  time  horizon 
(the  darker  the  grey,  the  more  commonly  it  occurred  on  the  installation). 
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Figure  11.  PCA  resultant  classes  that  occurred  within  Fort  Benning  in  the  2025  time  horizon. 


This  result  is  similar  to  the  Unsupervised  Classification  of  Figures  4  and  5 
that  show  categories  19  and  35  (respectively)  in  1990  migrating  to  Fort 
Benning  by  2025. 

The  most  common  class  on  Fort  Benning  for  years  2055  and  2085  does 
not  exist  in  the  1990  data  layer,  so  it  is  not  possible  to  show  where  the 
2055  and  2085  classes  come  from  in  the  1990  data.  It  is  worth  noting  that 
the  analysis  by  Unsupervised  Classification  also  did  not  generate  classes  in 
2055  and  2085  that  could  be  followed  in  the  1990  time  frame. 

4.6  The  Principally  Analogous  Multivariate  (PAM)  approach 

4.6.1  Introduction  to  the  PAM  approach 

Most  of  the  MVAs  previously  presented  are  based  on  the  classical  well 
documented  ISODATA  clustering  technique  from  the  Unsupervised 
Classification  method  (Ball  1965,  Jensen  1996,  Duda  2001).  In  carrying 
out  the  previous  approaches,  the  research  question  matured  as: 

If  a  similar  region  to  Fort  Benning  was  not  identified  in  the  time  period 
2055,  what  happened  to  it? 
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In  those  techniques  that  failed  to  identify  a  similar  patch  in  2055,  it  did 
not  seem  reasonable  to  conclude  that  the  patch  simply  vanished;  nor  did  it 
seem  reasonable  to  conclude  that  an  analogous  patch  did  not  exist.  It  must 
exist  somewhere,  but  some  of  the  other  MVA  techniques  were  restricted  by 
the  mathematical  procedure,  such  that  they  could  not  identify  what  the 
user  wanted  to  know.  The  Supervised  Classification  technique  did  produce 
usable  results,  apparently  because  the  user  in  that  situation  had  more 
control  over  the  process.  In  other  words,  the  user  provided  a  training 
sample  focused  on  identifying  characteristics  that  would  produce  given 
results,  to  then  see  if  those  results  would  better  fit  the  question.  Similarly, 
this  work  sought  a  more  direct  means  to  “let  the  system  know”  that  the 
results  of  interest  were  Fort  Benning  climate  migration. 

A  more  straightforward  technique  to  achieve  this  result  is  the  PAM 
approach.  This  technique  calculates  the  z-score  maps  for  each  of  the  19 
bioclimatic  concerns  (as  before).  Once  the  mean  of  an  area  of  interest  (in 
this  case  Fort  Benning)  is  determined,  this  mean  is  then  subtracted  from  the 
bioclimatic  map  of  interest.  The  smaller  the  resultant  value,  the  more 
similar  the  ensuing  areas  are  to  Fort  Benning.  If  you  add  up  all  the  resultant 
bioclimatic  maps,  then  the  smaller  the  output  value,  the  more  analogous 
areas  are  in  a  19-dimensional  bioclimatic  space.  Unlike  the  previous 
procedures,  this  technique  does  not  define  similar  locations;  instead  it 
defines  most  analogous  locations.  The  distinction  is  that  a  similar  location 
may  not  exist,  but  by  definition,  a  most  analogous  location  will  always  exist. 
The  advantage  of  this  procedure  is  that  it  ensures  an  analogous  location  in 
every  time  horizon.  PAM  makes  no  a  priori  assumptions  about  or 
requirements  on  the  input  data.  Within  the  results  are  also  contained  a 
rating  of  the  degree  of  analogousness  rather  than  a  simple  identification  of 
belonging  to  a  specific  class  as  in  the  other  techniques.* 

A  comparable  method  has  been  applied  to  ecosystem  migration  for 
military  installations  (Westervelt  2011).  In  that  analysis  the 
multidimensional  distance  between  the  target  signature  (representing  the 
future  state)  and  the  signature  at  each  location  in  the  “current”  data  was 


*  Some  software  packages  allow  another  map  output  (called  the  threshold  map)  that  does  rate  the 
degree  that  each  particular  cell  belongs  to  its  assigned  class. 
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calculated.  A  grey-scale-like  color  table  was  used  to  show  the  results.  The 
darker  areas  represented  more  similar  locations. 

4.6.2  Applying  the  PAM  approach  at  Fort  Benning 

When  applied  to  for  the  Fort  Benning  study  area  for  the  1990  time 
horizon,  this  procedure  results  in  a  19-dimensional  bioclimatic  map 
(Figure  12a).  The  darker  the  area,  the  more  analogous  that  location  is  to 
Fort  Benning.  Fortunately,  Fort  Benning  presents  as  the  darkest  location 
while  a  patch  to  the  east  is  also  comparable,  but  nowhere  else  in  1990. 
Applying  the  same  procedure  with  the  same  medians  to  the  similar  data 
for  2025  produces  slightly  different  results  (Figure  12b),  for  2055  (Figure 
12c)  and  for  2085  (Figure  I2d). 

These  analyses  show  that  the  current  bioclimatic  area  at  Benning  will 
slowly  migrate  about  210  km  to  the  north  east  over  the  next  90  years.  The 
objective  of  this  research  is  the  reverse  question: 

From  where  does  the  future  bioclimatic  character  of  Fort  Benning  come? 

The  same  technique,  modified  to  address  this  query,  generates  a  slightly 
different  set  of  images  (Figure  13). 

In  this  series  it  can  be  seen  that  the  bioclimatic  parameters  that  move  into 
Fort  Benning  come  from  a  region  about  56  km  southeast  of  the  installation 
in  2025,  then  about  93  km  south  in  2055  and  finally  138  km  south  in 
2085.  It  is  interesting  to  note  that  the  area  analogous  to  Fort  Benning  in 
1990  to  its  east  switches  sides  to  the  west  in  the  2055  and  2085  images.  In 
fact,  the  analogous  area  in  2085  in  central  Mississippi  is  just  as  strongly 
related  to  Fort  Benning  as  is  the  patch  to  the  south.  This  suggests  that 
there  are  possibly  two  bioclimatic  gradients:  one  marches  from  south  to 
north,  the  other  from  east  to  west.  The  east-to  west-march  also 
corresponds  with  the  results  shown  in  Figure  12,  which  indicate  that  Fort 
Benning’s  current  climate  will  relocate  toward  the  east.  Also,  invasive 
species  that  might  be  a  problem  to  Fort  Benning  are  equally  likely  to  come 
from  south-central  Mississippi  as  from  southwest  Georgia. 
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Figure  12.  The  19-dimensional  bioclimatic  PAM  analysis  for  the  Fort  Benning  study  area. 
These  analyses  represent  the  migration  of  the  climate  away  from  the  installation.  The  darker 
the  area,  the  more  analogous  that  location  is  to  Fort  Benning.  The  time  horizons  are: 
a  =  1990,  b  =  2025,  c  =  2055,  and  d  =  2085. 


c.  d. 

Figure  13.  Locations  from  “where  the  future  bioclimatic  character  of  Fort  Benning  come.”  The 
time  horizons  are:  a  =  1990,  b  =  2025,  c  =  2055,  and  d  =  2085. 
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The  Eufaula  National  Wildlife  Refuge  (Figure  14)  lies  about  50  km  to  the 
south  of  the  installation  and  falls  within  some  of  the  patches  that  migrate 
north  to  Fort  Benning  as  just  discussed.  By  2025  the  climate  around  Fort 
Benning  will  be  tending  to  look  like  what  is  currently  within  the  refuge.  A 
healthy,  growing  population  of  threatened  American  alligators  exists  in 
refuge  wetlands  (USFWS  Undated).  Resident  wildlife,  including  bobwhite 
quail,  wild  turkey,  beaver,  fox,  raccoon,  opossum,  bobcat,  swamp  and 
cottontail  rabbits,  and  white-tailed  deer  are  present  in  high  numbers. 

The  northerly  migration  of  the  bioclimatic  concerns  suggest  that  in  the 
next  period,  about  2115  (for  which  data  were  unavailable),  Fort  Benning’s 
climate  will  mimic  a  location  in  the  Florida  Panhandle  that  is  occasionally 
called  the  Gulf  Coast  Sandhills.  At  that  time  period,  not  only  will  the 
climate  be  similar,  but  the  soil  character  will  also  be  similar.  Therefore 
about  the  year  2115,  invasives  from  the  Florida  Panhandle  will  be  highly 
capable  of  taking  up  their  residence  at  Fort  Benning. 

Figure  14.  Upland  landscape  near  the  entrance  to  the  Eufaula  National  Wildlife  Refuge. 
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The  data  in  Table  6  outline  the  characteristics  of  these  regions  for  each  of 
the  time  periods.  Table  6  lists  the  values  and  amount  of  change  that  can  be 
extracted  from  the  areas  defined  above  for  each  of  the  19  bioclimatic 
concerns.  This  data  answers  the  research  question  posed  at  the  outset  of 
the  report,  i.e.,  What  will  the  climatic  factors  in  the  future  be  at  the 
installation?”  A  similar  table  could  be  generated  from  the  Supervised 
Classification  technique,  but  not  the  Unsupervised  Classification  or 
Principal  Components  because  they  did  not  generate  results  that  could  be 
followed  in  2055  and  2085. 

4.6.3  Migration  of  Fort  Benning  forest  types  —  lowlands 

During  the  late  1990s,  the  Fort  Benning  Land  Management  Branch 
compiled  a  digital  Forest  Inventory  map  of  the  forest  types  at  the 
installation.  This  map  provides  a  snapshot  of  the  woodland  distributions 
at  roughly  the  time  (the  late  1990s)  as  the  base  climate  data  (1990).  It 
would  be  useful  if  it  were  possible  to  use  this  map  and  a  method  similar  to 
the  PAM  analysis  in  the  last  section  to  follow  the  location  of  incoming 
bioclimatic  changes  to  Fort  Benning.  The  map  is  used  to  find  the 
bioclimatic  parameters  for  just  the  tree  types  rather  than  the  entire 
installation.  This  set  of  parameters  is  then  applied  to  the  temporal  data  to 
follow  changes  over  time.  First  the  inventory  was  used  to  generalize  the 
species  into  less  than  a  dozen  major  types  (Figure  15).  That  inventory  was 
narrowed  into  just  those  wetland  types  (blue)  within  a  circle  centered  on 
the  Lower  Upatoi  Creek  that  crosses  Fort  Benning’s  northwest  quarter 
(also  Figure  15).  The  result  (Figure  16)  shows  the  1990  bioclimatic 
distribution  of  cells  that  coordinate  closely  with  the  lowland  tree  species 
similar  to  the  Lower  Upatoi  Creek. 

Some  sensitivity  tests  were  carried  out  on  the  results.  First,  the  median 
values  were  generated  at  the  resolution  of  the  forest  inventory 
(0:00:00.323629  degree  or  10  meters)  and  the  map  generated.  Then  the 
median  values  were  developed  at  the  bioclimatic  data  resolution  (0:00:30 
degree  or  about  925  meters)  and  the  map  was  generated. 


Table  6.  Characteristics  of  the  PAM  derived  regions  for  each  of  the  time  periods  for  each  of  the  19  bioclimatic  concerns 
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Figure  15.  Generalize  species  from  the  Fort  Benning  Forest  Inventory  map.  Wetland  types  are 
in  blue.  The  target  areas  are  the  blue  within  the  circle  centered  on  the  Lower  Upatoi  Creek. 


Figure  16.  The  1990  bioclimatic  distribution  of  cells  that  coordinate  closely  with  the  lowland 
tree  species  similar  to  the  Lower  Upatoi  Creek  from  Fig.  25. 


ERDC/CERL  TR-15-2 


37 


Although  the  median  values  varied  slightly,  the  resultant  cells  chosen  were 
almost  the  same.  In  another  set  of  tests,  a  small  constant  was  subtracted 
from  all  19  median  values.  (Mean  values  all  varied  slightly  about  the  value 
127.)  There  was  no  change  if  the  constant  subtracted  was  0.001,  but  an 
entirely  different  result  was  obtained  when  the  constant  subtracted  reached 
0.005.  Thus  the  result  is  not  highly  sensitive  to  changes  of  resolution,  but 
can  be  greatly  altered  by  changing  of  the  mean  values  by  as  little  as  0.004%. 

The  PAM  analysis  was  applied  to  the  different  time  horizons.  Figure  17 
shows  the  resultant  situation  for  2025.  The  darker  grey  areas  show  all  of 
the  distribution  previously  chosen  for  all  of  the  installation  (as  shown  in 
Figure  13b)  while  the  lighter  grey  areas  show  just  the  lowlands  chosen  in 
this  analysis.  The  lighter  areas  clearly  tend  to  congregate  about  lowlands, 
either  along  the  Chattahoochee  River  or  particularly  along  the  ravines  in 
the  Lumpkin  region  about  50  km  to  the  south  of  Fort  Benning. 

Next  the  analysis  was  constrained  by  integrating  the  geomorphology.  The 
enhanced  algorithm  penalized  upland  locations  (using  a  geomorphology 
dataset  generated  from  a  topography  data  layer  with  a  higher  resolution 
than  the  bioclimatic  data).  Figure  17b  integrates  the  geomorphology  to 
refine  the  analysis.  In  this  case,  the  algorithm  provides  a  bias  to  areas  that 
are  bioclimatically  appropriate  and  also  lowlands.  Henceforth  in  this 
section,  all  graphics  will  show  the  results  with  the  geomorphology 
integrated  into  the  image. 

By  2055  (Figure  18),  all  of  the  lowland  habitat  on  the  Chattahoochee  west 
of  Fort  Benning  has  disappeared,  the  prime  Lumpkin  locations  50  km  to 
the  south  remain  and  the  areas  120  km  to  the  southeast  near  Albany  are 
much  more  prominent.  The  two  major  locations  from  which  the  lowland 
vegetation  will  come  (the  Chattahoochee  and  the  upper  Hodchodkee  and 
Pataula  Creek  areas)  are  largely  the  same  in  2025  as  in  2055  with  the 
caveat  that  they  move  south  by  a  few  kilometers  in  both  cases. 

In  about  2085,  the  prime  areas  will  have  migrated  even  further  south,  to 
near  the  Florida  state  line  (Figure  19).  Although  these  areas  still  tend  to 
follow  the  lower  neighborhood,  at  this  time  period,  some  inter-ravine 
areas  are  also  included.  This  region  represents  a  topographic  character 
that  is  less  well  defined;  therefore  the  landform  constraints  are  less  clear. 
The  most  prominent  focus  is  no  longer  near  the  Chattahoochee  River,  but 
along  the  Flint  River  and  Spring  Creek  to  the  east  of  the  Chattahoochee. 
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Figure  17.  PAM  analysis  for  the  distribution  of  lowland  tree  vegetation  change  in  2025.  Right: 
The  darker  grey  areas  show  all  of  distribution  previously  chosen  for  all  of  the  installation  (as 
shown  in  Fig.  13,  upper  right)  while  the  lighter  grey  areas  show  just  the  lowlands  chosen  in 
this  analysis.  (All  further  graphics  in  this  section  will  use  this  convention.)  Left:  Same  analysis 
with  geomorphology  integrated  into  the  result. 
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Figure  19.  Lowland  climate  source  for  Fort  Benning  in  2085. 


As  mentioned  earlier,  there  is  an  interesting  additional  locus  of  similar 
climatic  habitat  in  central  Alabama  south  of  Selma.  Focusing  in  on  this 
area  reveals  a  similar  type  fit  between  the  chosen  lowland  areas  (Figure 
20)  and  river  locations.  This  locus  however  is  even  brighter  (meaning  the 
fit  is  better)  and  more  compact  than  the  previous  one  near  the  Florida 
boundary.  Furthermore,  it  is  similarly  poorly  constrained  by  topography. 
This  central  Alabama  location  is  currently  highly  characteristic  of  the 
climate  in  the  Lower  Upatoi  River  that  Fort  Benning  can  expect  to 
experience  in  2085. 

4.6.4  Migration  of  Fort  Benning  forest  types  -  uplands 

This  section  investigates  whether  it  is  possible  to  distinguish  a  different 
type  of  vegetation  migration  over  time.  The  type  of  vegetation  chosen  was 
upland  pine  types  (shown  as  yellow  in  Figure  21)  within  a  circle  centered 
on  the  southeast  corner  of  Fort  Benning  (also  shown  in  Figure  21).  Figure 
22  shows  the  1990  bioclimatic  distribution  of  cells  that  coordinate  closely 
with  the  upland  pine  tree  species  in  the  southeast  corner  of  Fort  Benning. 
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Figure  20.  Alternative  lowland  climate  source  from  Alabama  for  Fort  Benning  in  2085. 


Figure  21.  Generalize  species  from  the  Fort  Benning  Forest  Inventory  map.  Upland  types  are 
in  yellow.  The  target  areas  are  the  yellow  within  the  circle  centered  on  the  southeast  corner  of 

the  installation. 
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Figure  22.  The  1990  bioclimatic  distribution  of  cells  that  coordinate  closely  with  the  upland 
tree  species  similar  to  those  in  the  southeast  corner  of  Fig.  31. 


Figure  23  shows  the  resulting  upland  pine  distribution  for  2025.  The 
darker  grey  areas  show  the  full  extent  of  the  2025  distribution  (as 
previously  shown  in  Figure  13  upper  right)  while  the  lighter  grey  areas 
show  just  the  pixels  chosen  in  this  analysis.  This  time,  the  lighter  areas 
tend  to  be  congregated  along  the  uplands  between  the  ravines  near  the 
Chattahoochee  River,  particularly  along  the  uplands  between  the  ravines 
in  the  region  south  of  Lumpkin  about  50  km  to  the  south  of  Fort  Benning. 

By  2055,  the  climate  has  encouraged  the  upland  species  to  be  derived  from 
an  area  further  south  (Figure  24).  Although  the  uplands  near  Hodchodkee 
and  Pataula  Creeks  are  still  present  in  this  analysis,  it  is  their  lower  reaches 
(the  more  southerly)  that  are  now  favored.  Furthermore,  the  southeast  area, 
which  only  had  a  speckling  of  upland  representatives  in  2025,  is  now  much 
more  highly  represented  in  2055.  Note  that,  as  one  progresses  to  the  south, 
the  terrain  changes  from  the  rolling  hills  characteristic  of  Fort  Benning  to 
the  flatter  coastal  plains  where  the  drainage  is  more  poorly  defined  (Figure 
24).  Because  of  this  change  in  landform  characteristics,  the  corresponding 
distinction  between  uplands  and  lowlands  is  less  well  defined  and  this  is 
reflected  in  the  resulting  speckled  patterns  seen  in  this  Coastal  plains  area. 
Figure  25  shows  the  same  results  draped  over  the  topography. 
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Figure  23.  Fort  Benning  source  for  upland  pine  distribution  for  2025. 


Figure  24.  Fort  Benning  source  for  upland  pine  distribution  for  2055. 
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Figure  25.  Upland  pine  distribution  for  2055  (from  Fig.  24)  draped  over  the  topography  map. 


By  2085,  the  area  from  which  upland  vegetation  is  migrating  is  near  the 
Florida  border,  but  because  of  the  indeterminate  topography,  the  locations 
are  scattered  as  shown  in  Figure  26. 


4.6.5  Analysis  of  the  confusion  between  upland  and  lowland  plant 
locations 

It  is  apparent  in  Figure  26  that  the  distinction  between  the  upland  and 
lowland  species  is  less  clear.  Therefore,  a  research  question  of  interest  is: 

How  distinct  is  the  delineation  between  upland  and  lowland  sources  at  a 
specific  time  horizon?” 


A  sampling  technique  was  developed  to  determine  exactly  how  well  this 
analysis  separated  the  different  areas  identified  in  the  previous  two 
sections.  Figure  27  shows,  for  the  year  2025,  the  distribution  of  upland 
species  (colored  blue),  lowland  species  (green)  and  the  location  where 
both  upland  and  lowland  species  were  identified  in  the  same  cell  location 
(red).  Red  cells  show  the  amount  of  “confusion”  in  this  analysis.  There  are 
only  a  few  red  cells  along  the  Chattahoochee  River  near  Fort  Benning. 
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Figure  26.  Fort  Benning  source  for  upland  pine  distribution  for  2085. 


Figure  27.  The  distribution  for  the  year  2025  of  upland  species  (colored  blue),  lowland 
species  (green)  and  where  both  upland  and  lowland  species  were  identified  in  the  same  cell 

location  (red). 
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Table  7.  Upland  vs.  lowland  confusion  table  in  each  of  the  time  horizons 


Time  Horizon 

Only  Upland 

Only  Lowland 

Both  Upland 
and  Lowland 

Total  Number 
of  Cells 
Identified 

Percent 

Confused 

Cells 

2025 

26596 

31632 

527 

58755 

0.90% 

2055 

30826 

18773 

0 

49599 

0.00% 

2085 

48726 

14052 

0 

62778 

0.00% 

A  similar  analysis  was  carried  out  for  the  other  two  time  horizons,  2055 
and  2085.  Table  7  lists  the  results  of  all  three  analyses;  the  data  show  that 
the  distributions  are  almost  mutually  exclusive.  (The  worst  case  is  less 
than  1%  confusion  for  2025.)  Since  the  visual  distribution  of  forest  types 
for  the  2085  distribution  shown  in  Figure  24  was  considered  less  distinct, 
the  data  in  Table  7  indicate  that  the  method  works  extremely  well.  This 
analysis  is  good  evidence  that  sub-installation  distributions  of  vegetative 
types  can  be  reliably  identified  as  climate  changes  occur  over  time. 

4.6.6  Invasive  species  based  on  the  Primarily  Analogous  Multivariate 
approach 

Invasive  Species  (IS)  can  come  from  two  sources: 

•  nearby  similar  areas 

•  distant  locations  with  similar  characteristics. 

After  exploring  the  potential  sources  of  IS  from  nearby  areas,  one  may  also 
consider  where  analogous,  distant  locations  exist.  Assuming  that  species 
from  distant  locations  that  have  a  highly  similar  environment  are  most 
likely  to  successfully  establish  themselves  around  Fort  Benning,  the 
research  question  can  now  be  rephrased  as: 

What  are  the  bioclimatic  regions  elsewhere  in  the  world  that  are  most 
analogous  to  the  region  around  Fort  Benning? 

Since  climate  change  will  alter  the  current  most  analogous  area  to  a 
vicinity  that  might  even  have  a  different  ecosystem,  it  is  necessary  to  look 
both  at  analogous  areas  in  the  current  time  horizon  as  well  as  those  in  the 
future.  The  technique  used  here  is  comparable  to  the  regional  Primarily 
Analogous  Multivariate  approach,  but  applied  globally.  Figure  28  shows 
that,  after  the  analysis,  Fort  Benning  is  fortunately  similar  to  itself  (very 
dark  areas  in  the  figure)  and  to  the  general  region  of  the  southeastern 
United  States.  There  is  no  other  extensive  region  in  the  world  that  is  so 
akin  to  the  southeastern  United  States. 
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Figure  28.  The  resultant  PAM  bioclimatic  regions  elsewhere  in  the  world  that  are  most  similar 
in  1990  to  the  region  around  Fort  Benning  (the  darker,  the  more  similar). 


The  hardiest  invasives  at  Fort  Benning  are  expected  to  be  from  the 
southeastern  United  States.  However,  some  more  distant  regions  are 
analogous.  Specifically,  in  order  of  decreasing  coverage,  the  similar  areas 
are  in: 

•  South  America  along  the  border  of  Uruguay  and  Argentina  (the  Parana 
River  Valley) 

•  Along  the  eastern  edge  of  Australia 

•  Southern  Europe  and  northern  Mediterranean 

•  East  central  China. 

Although  the  Parana  River  area  is  more  extensive,  the  volume  of  trade 
from  Europe  makes  this  a  more  likely  source  of  foreign  species.  South 
American  species  may  have  a  difficult  task  adapting  to  North  America 
because  the  seasons  are  reversed.  This  is  not  true  of  European  or  Asian 
species.  These  are  the  areas  from  which  invasives  are  most  likely  to  come 
that  end  up  in  the  southeast  United  States. 

The  companion  analysis  for  2085  is  instructive  (Figure  29).  It  shows 
shrinkage  of  similar  areas  from  Europe.  On  the  other  hand,  the  region 
around  the  Parana  has  increased  as  has  the  region  in  China.  Small  new 
areas  can  now  be  seen  intensifying  in  Northern  Pakistan.  Australia  has 
become  less  intensely  similar,  but  that  which  exists  has  spread  over  much 
more  of  the  continent.  Finally  the  Sahel  region  of  Africa  has  now  emerged 
prominently  whereas  it  was  barely  represented  in  1990. 
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Figure  29.  The  PAM  analysis  for  2085  for  the  bioclimatic  regions  that  are  most  similar  to  the 

region  around  Fort  Benning. 


This  suggests  that,  with  increased  trade,  invasives  from  the  Parana  basin 
will  become  more  significant  to  Fort  Benning.  China,  Australia,  the  Sahel 
and  Northern  Pakistan  will  provide  more  adopted  IS  to  the  installation 
than  they  do  now.  However,  invasives  from  the  southeast  region  of  the 
United  States  are  still  likely  to  have  an  easier  time  finding  a  foothold  at 
Fort  Benning. 

The  US  Army  Floristic  Inventories  database  (HQUSACE  2007)  lists  101 
exotic  vascular  plant  species  at  Fort  Benning.  Of  those,  the  species  that  are 
most  often  the  focus  of  attention  include:  Perennial  Kudzu  ( Pueraria 
lobata ),  Chinese  Bush  Clover  ( Lespedeza  cuneata ),  Oriental  Bittersweet 
( Celastrus  cn'biculatus),  Butternut  Canker  ( Sirococcus  clavigignenti- 
juglandacearum )  and  Cogon  grass  (Imp  er  at  a  cylindrica).*  Kudzu  is 
native  to  the  Indian  subcontinent;  Chinese  Bush  Clover  is  from  Asia  and 
eastern  Australia;  Oriental  Bittersweet  is  native  to  East  Asia;  Butternut 
Canker  is  possibly  from  Asia;  and  Cogon  grass  is  native  to  east  and 
southeast  Asia,  India,  Micronesia,  Australia,  and  eastern  and  southern 
Africa.  One  of  the  listed  exotics  is  South  American  mock  vervain 
0 Glandularia  pulchella ),  which  is  native  to  Brazil,  Argentina,  and 
Uruguay. 

These  similar  areas  are  all  represented  in  Figures  28  and  29.  The 
prevalence  of  Asian  invasives  suggests  strong  historical  trade  connections. 
The  lack  of  problem  species  from  Europe  is  interesting  to  note.  It  has  been 


*  http://workshops.serdp- 

estcp.org/serps/docs/Southeastern_RegionaLConservation_Research_Projects.pdf 
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suggested  that  European  invaders  have  been  in  the  United  States  long 
enough  that  they  are  now  considered  naturalized.  No  references  were 
found  to  support  this  suggestion. 

Table  8  lists  the  proposed  exotic  pest  plant  species  for  the  State  of  Georgia 
(Murphy  2005).  All  of  them  come  from  Asia,  usually  China. 

The  analysis  in  this  section  suggests  that  more  invasives  may  be  expected 
to  arrive  in  the  southeastern  United  States  from  China  due  to  climate 
change. 


Table  8.  Proposed  exotic  pest  plant  species  for  Georgia. 


Common  Name 

Scientific  Name 

Native  Land* 

Kudzu 

Pueraria  montana  (Lour.)  Merr. 

Indian  Subcontinent 

Chinese  privet 

Ligustrum  sinense  Lour. 

China,  Taiwan  and  Vietnam 

apanese  honeysuckle 

Lonicera  japonica  Thunb. 

China,  Japan  and  Korea 

Hydrilla 

Hydrilla  verticillata  (L.  f.)  Royle 

Asia,  Europe,  Africa  and 

Australia 

Chinese  tallow  tree 

Triadica  sebifera  (L.)  Small 

Eastern  China,  Taiwan,  and 
Japan 

Nepalese  browntop 

Microstegium  vimineum  (Trin.)  A. 
Camus 

Iran  in  the  west,  east  to  China 

Bamboo 

Phyllostachys  aurea  Carr,  ex  A.& 

C.  Riviere 

China 

Autumn  olive  (silverberry) 

Elaeagnus  umbellata  Thunb. 

Eastern  Asia  to  Japan 

Chinese  wisteria 

Wisteria  sinensis  (Sims)  DC. 

China 

Mimosa 

Albizia  julibrissin  Durazz. 

Southwestern  to  eastern  Asia 

'Native  Lands  column  is  compiled  from  Wikipedia  entry  for  each  species. 
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5  Comparison  of  the  Different  Multivariate 
Approaches 

5.1  The  evaluation  of  the  different  techniques 

This  research  investigated  the  application  success  and  versatility  of  six 
different  multivariate  approaches.  Two  approaches,  (CCA)  and  an  SMAP 
estimation,  provided  no  usable  signature  files  after  expending  a  great  deal 
of  effort.  Table  9  lists  a  summary  of  the  evaluation  of  the  six  MVAs  applied 
in  this  work. 


Table  9.  Evaluation  of  six  MVAs. 


Approach 

Pros 

Cons 

Unsupervised 

Classification 

•  Well  accepted  and  understood  procedure 

•  Requires  minimal  data  preparation 

•  Requires  minimal  additional  user  input 

•  Requires  few  assumptions  on  the  part  of 
the  user 

•  Results  are  roughly  consistent 
independent  of  the  number  of  layers 
input 

•  Results  are  not  tied  to  known  attributes 
on  the  land 

•  Lack  of  output  control 

•  Did  not  result  in  patches  that  can  be 
followed  through  time 

•  User  cannot  control  well  the  character  of 
the  output  -  this  is  mostly  determined  by 
the  algorithm 

Supervised 

Classification 

•  Well  accepted  and  understood  procedure 

•  Requires  training  map  -  you  know  the 
meaning  of  the  resultant  areas 

•  Resulted  in  patches  that  can  be  followed 
through  time 

•  User  cannot  control  well  the  character  of 
the  output  -  this  is  mostly  determined  by 
the  algorithm 

•  Requires  training  map 

Principal 

Components 

•  Sophisticated  analysis  approach 

•  Reforms  the  data  to  make  them  more 
acceptable  to  classifying  programs 

•  Level  of  variance  in  data  numerically 
assured 

•  Lack  of  output  control 

•  Did  not  result  in  patches  that  can  be 
followed  through  time 

•  Data  are  tertiary  in  generation  -  there  is 
a  chance  the  original  information  value 
has  been  lost 

Primarily 

Analogous 

Multivariate 

Approach 

•  Easy  to  conceptually  describe 

•  Easy  to  carry  out 

•  It  puts  no  restrictions  on  the  input  data 

•  Guaranteed  to  result  in  patches  that  can 
be  followed  through  time 

•  Versatile  -  allows  sub-installation  and 
global  level  analyses 

•  Not  sophisticated  in  concept 

Canonical 

Components 

Analysis 

N/A 

Produced  no  results 

SMAP  Estimation 

N/A 

Produced  no  results 
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5.2  The  effect  of  the  Benning  Hole  on  the  results 

Chapter  3  discussed  an  anomaly  in  the  data,  called  the  “Benning  Hole,”  in 
which  “patches  50  kilometers  to  the  east  of  the  installation  are  expected  to 
be  more  similar  to  the  installation  than  areas  on  it  or  more  immediately 
adjacent.”  Evidence  suggests  that  this  potential  anomaly  influenced  the 
results  of  these  analyses. 

Several  of  the  outputs  show  a  clear  east-west  bias  (as  shown  in  Figure  8  in 
particular,  but  also  in  Figures  4,  9, 10,  24,  and  25).  In  other  words,  several 
of  the  analyses  support  the  hypothesis  that  an  east-west  component  seems 
to  exist.  More  specifically,  the  analysis  indicates  that  “the  current 
Bioclimatic  at  Benning  will  go”  on  a  West  to  East  track  (Figure  12)  instead 
of  a  south  to  north  track.  On  the  other  hand,  the  elevation  rises  quickly  to 
the  north  of  the  installation  that  a  purely  northerly  migration  would  not  be 
expected  on  that  basis  alone.  There  is  evidence  (but  no  clear  proof)  for  the 
influence  of  the  hole  in  the  results. 

On  the  other  hand,  most  of  analyses  show  the  new  Fort  Benning  climate 
migrating  from  south  to  north  (particularly  Figure  13),  roughly 
perpendicular  to  the  trend  of  the  values  associated  with  the  hole.  The  final 
conclusion  on  the  issue  is  that,  if  the  hole  influences  the  result,  it  is  not 
great  enough  to  overpower  the  overall  migration  of  the  bioclimate  in  a 
direction  perpendicular  to  that  which  the  hole  should  produce. 
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6  Conclusions  and  Recommendations 

6.1  Conclusions 

Changes  in  climatic  parameters  are  important  in  that  they  affect  the 
military’s  ability  to  perform  its  national  defense  mission  to  manage 
military  lands  within  the  continental  United  States  (CONUS).  The  military 
must  have  the  ability  to  predict  climatic  changes  on  specific  installations. 
To  support  this  need,  the  research  has  outlined  and  evaluated  the 
application  of  MVA  techniques  to  predict  climatic  changes  on  a  specific 
Army  installation,  Fort  Benning  GA. 

To  accomplish  the  evaluation,  a  database  was  created  of  19  bioclimatic 
parameters  that  represented  a  consensus  of  six  GCMs  and  three  future 
scenarios  for  each  of  the  time  periods  centered  around  the  dates  1990, 
2025,  2055,  and  2085.  Once  in  place,  this  data  helped  to  form  a  solid  base 
that  made  it  possible  to  reliably  track  the  predicted  changes.  Several  MVA 
techniques  were  applied  to  the  data  and  the  outcomes  were  compared. 

The  results  of  almost  all  the  investigations  indicated  that  predicted  climate 
change  was  found  to  be  great  enough  so  that  almost  no  place  over  the 
entire  installation  would  be  the  same  in  2025  as  it  had  been  in  1990. 

Of  the  six  MVA  techniques  evaluated,  four  provided  usable  results: 

•  Unsupervised  Classification 

•  Supervised  Classification 

•  Principal  Components 

•  Primarily  Analogous  Multivariate. 

The  Supervised  Classification  and  the  Primarily  Analogous  approaches 
produced  usable  predictions  for  all  four  time  periods.  The  Unsupervised 
Classification  and  Principal  Components  techniques  produced  usable 
predictions  for  only  2025. 

The  evaluation  found  that  the  Primarily  Analogous  Multivariate  approach 
developed  during  this  research  clearly  distinguished  itself  from  the  other 
five  approaches  in  that  it: 

•  was  easy  to  conceptualize 

•  was  easy  to  carry  out 

•  required  no  restrictions  on  the  input  data 
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•  guaranteed  that  results  could  be  generated  at  each  time  horizon 

•  was  highly  versatile  in  that  it  successfully  supported  application  from 
the  sub-installation  to  the  global  scales. 

Of  the  six  approaches  evaluated,  two  (the  Supervised  Classification  and 
PAM)  could  fulfill  the  primary  objective  of  this  research,  to  determine  “the 
climatic  factors  in  the  future  at  the  installation.” 

Of  the  six  approaches  evaluated,  one  (PAM)  could  fulfill  the  secondary 
objective  of  this  research,  to  determine  whether  “an  analysis  at  the  sub¬ 
installation  level  can  be  successfully  accomplished  with  the  currently 
available  climatic  data.” 

6.2  Recommendations 

As  a  result  of  this  investigation,  it  is  recommended  that  research  continue 
to  develop  the  capability  to  anticipate  and  predict  how  the  bioclimate  will 
change  over  time  at  a  particular  installation.  Furthermore,  it  is 
recommended  to  promote  the  application  of  the  Primarily  Analogous 
Multivariate  approach  to  analyze  climate  change  at  the  sub-installation 
level.  This  simplest  of  the  six  evaluated  approaches  always  produced 
usable,  reliable,  detailed  outputs,  and  exhibited  the  greatest  versatility  in 
being  applied  to  different  issues  at  scales  from  the  sub-installation  to  the 
global  level. 
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Acronyms 

Term 

CCA 

CEQ 

CERL 

Cl  AT 

CONUS 

DoD 

ERDC 

ERDC-CERL 

GCM 

GFDL 

GISS 

IEEE 

INRMP 

IPCC 

IS 

MVA 

NASA 

NCAR 

NEPA 

NOAA 

NSN 

OMB 

PAM 

PCA 

QDR 

SMAP 

SR 

SRES 

TAR 

TR 

URL 

US 

USA 

WWW 


and  Abbreviations 


Definition 

Canonical  Components  Analysis 
Council  on  Environmental  Quality 
Construction  Engineering  Research  Laboratory 
International  Centre  for  Tropical  Agriculture 
Continental  United  States 
US  Department  of  Defense 

US  Army  Engineer  Research  and  Development  Center 

Engineer  Research  and  Development  Center,  Construction  Engineering 
Research  Laboratory 

Global  Climate  Model 

Geophysical  Fluid  Dynamics  Laboratory 

Goddard  Institute  for  Space  Studies 

Institute  of  Electrical  and  Electronics  Engineers 

Integrated  Natural  Resources  Management  Plans 

Intergovernmental  Panel  on  Climate  Change 

Invasive  Species 

Multivariate  Analysis 

National  Aeronautics  and  Space  Administration 

National  Centre  for  Atmospheric  Research 

National  Environmental  Policy  Act 

National  Oceanic  and  Atmospheric  Administration 

National  Supply  Number 

Office  of  Management  and  Budget 

Primarily  Analogous  Model 

Principal  Components  Analysis 

Quadrennial  Defense  Review 

Sequential  Maximum  A  Posteriori 

Special  Report 

The  Special  Report  on  Emissions  Scenarios 

Third  Assessment  Report 

Technical  Report 

Universal  Resource  Locator 

United  States 

United  States  of  America 

World  Wide  Web 
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